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We report first-principles calculations of acoustic surface plasmons on the (0001) surface of Be, as obtained 
in the random-phase approximation of many-body theory. The energy dispersion of these collective exci- 
tations has been obtained along two symmetry directions. Our results show a considerable anisotropy of 
acoustic surface plasmons, and underline the capability of experimental measurements of these plasmons to 
map the electron-hole excitation spectrum of the quasi two-dimensional Shockley surface state band that is 
present on the Be(OOOl) surface. 
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1 Introduction 

The existence of metal surfaces is well known to yield the appearance of conventional surface collective 
electronic excitations, surface plasmons, predicted by Ritchie 50 years ago H]. Thereafter, surface plas- 
mons were widely investigated both experimentally and theoretically 12] [3] |4l |5] |6l 13 • Another kind of 
collective excitations that are present at metal surfaces are the so-called multipole surface plasmons pre- 
dicted theoretically fSl and later detected experimentally 19] [TO]. The energy of both these plasmons is a 
fraction of the bulk plasmon energy, ujp, typically (0.7 — 0.8)wp for small momenta. 

A qualitatively different surface collective excitation, called acoustic surface plasmon, has recently 
been predicted theoretically 111] and observed experimentally using angle-resolved electron energy loss 
spectroscopy |12J. Acoustic surface plasmons are low-energy collective excitations that are confined to 
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solid surfaces where a partially occupied quasi-two dimensional (2D) surface-state band coexists with 
the underlying three-dimensional (3D) continuum, as occurs in the case of Be(OOOl) |11 1 and the (111) 
surfaces of the noble metals Cu, Ag, and Au [,13]. They are called acoustic because their energy dispersion 
exhibits sound-like linear behaviour at low momenta, contrary to the square-root dispersion expected for a 
pure 2D electron gas [T?]. 

The calculations that led to the prediction of acoustic surface plasmons ifTTl assumed translational in- 
variance in the plane of the surface and incorporated the electron dynamics in the direction perpendicular 
to the surface through the use of a one-dimensional (ID) model potential that describes the main features 
of the surface band structure LISJ . Although acoustic plasmons had previously been only expected to exist 
for spatially separated plasmas lfT6l[T7l . ID model potential self-consistent calculations demonstrated that, 
contrary to expectations, acoustic plasmons should indeed be present at bare metal surfaces where a 2D 
surface-state band coexists with the underlying 3D continuum. 

In this paper, we report first-principles calculations of acoustic surface plasmons on the (0001) surface 
of Be, as obtained in the random-phase approximation of many-body theory. First-principles calculations 
along the TK symmetry direction were reported before |„12J. but here we extend those calculations to 
other symmetry directions and focus on the discussion of the anisotropy of acoustic surface plasmons. 
Our calculations show that the anisotropy is considerable and underline the capability of experimental 
measurements of these plasmons to map the electron-hole excitation spectrum of the quasi two-dimensional 
Shockley surface state band that is present on the Be(OOOl) surface. 



2 Calculation details 

The rate at which a frequency-dependent external potential (j)™*{r, uj) generates electronic excitations in a 
many-electron system can be obtained, within lowest-order perturbation theory, as follows 161 

u;H = -2Im y"drni„d(r,u;)r"*(r,w), (1) 

where nind(r, i^) represents the induced electron density. In the case of a periodic solid surface, one can 
write Q 

SBZ 

w{lo) — ^ w{q, uj), (2) 
q 

where the sum extends over the surface Brillouin zone (SBZ) and w(q, ui) denotes the rate at which the 
external potential generates electronic excitations of frequency uj and parallel wave vector q: 



w;(q,w) = --^Im J dz nind(2, q, w) q, w), 



(3) 



nind{z, q, and cf)'^^* {z,q_, lo) representing 2D Fourier transforms of the induced electron density nind(i", 
and the external potential (^"^^'(r, lo), respectively, and A being the normalization area. If the 2D Fourier 
transform of the external potential is of the form 



'(0,q,cj) = — 



27r 



(4) 



the rate w(q, oj) can be expressed in the following form: 



w(q,w) 



47r _ 
— Im5(q,w), 
qA 



(5) 
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where g(q, uj) represents the so-called surface response function ifTSl 

2tt f 

gi'^,^) = — z- J dzdz'xG=oG'=o(^'^''^'^)' (6) 

Xg=og'=o('^' -^ 'q'^) being the Fourier components of the density response function of the interacting 
many-electron system. In the random-phase approximation of many -body theory, these Fourier coefficients 
are easily obtained from the Fourier coefficients X^-^' {^i ^'7 ^) of the noninteracting density response 
function by solving the following matrix equation: 

XGG' (q, w) = xgg' (q> ^) + ^ xggi (q> ^) '^d (q) xgig- (q, w) , (7) 

Gi 

where vg (q) ~ \ q+G[^ represent the Fourier coefficients of the bare Coulomb electron-electron interac- 
tion. 

For positive frequencies, the imaginary part of the Fourier coefficients Xgg' (q^ ^) of '^^e noninteracting 
density response function is easily obtained from the spectral function (q, uj), as follows 

im[xGG'(q7^)] = -^5^G'(q7^)7 (8) 

where 

^ SBZ occ unocc 

5GG'(q,-) = ^EEE(^kJ-~'*^'^°^^iW) 



X (V'k+q„'|e'('^'+''''"^IVl^„)'5(ek„-£k+q,,+c-), (9) 

the factor 2 accounting for spin. The sum over n in') runs over occupied (unoccupied) states of the sur- 
face band structure for 2D wave vectors k in the SBZ. We take and tp^^^ to be the eigenvalues and 
eigenfunctions of a first-principles single-particle Kohn-Sham (KS) Hamiltonian of density-functional the- 
ory (DFT). For the evaluation of the real part of the Fourier coefficients Xgg' (q' ^) of noninteracting 
density response function, we perform a Hilbert transform of the corresponding imaginary part. 

For a description of the (0001) face of hep Be from first principles, we employ a repeated-slabs geom- 
etry. The slabs periodically repeat in the direction normal to the (0001) surface and every slab consists 
from 24 atomic Be layers. The distance between slabs is chosen to correspond to 8 Be atomic layers. For 
the evaluation of the eigenvalues and eigenfunctions entering Eq. (|9]l we solve the Kohn-Sham equation of 
DFT self-consistently with the use of norm-conserving pseudopotentials ||T9l . the local-density approxima- 
tion (LDA) for exchange and correlation, and the Perdew-Zunger parametrization of the Ceperley-Alder 
diffusion Monte Carlo energies of a uniform electron gas ||20l |2TI . We expand the single-particle wave 
functions i/^ ^ (r) in a plane wave basis, with reciprocal vectors G up to an energy cut-off of 20 Ry. For 
the lattice constants, we take the experimental values a — 2.285 A and c = 3.585 A |22|. The experimental 
value ll23l for the relaxation of surface atomic layers has been taken into account. The sum over 2D wave 
vectors k in Eq. (|9|l is taken to run over a 108 x 108 mesh, which corresponds to 1 1664 points in the SBZ. 
The energy band summation includes all bands up to an energy cut-off of 50 eV above the Fermi level. The 
integrations involved in the evaluation of the surface response function of Eq. (|6]l are performed from — c/2 
up to c/2, where c is the unit cell size in the z-direction. Finally, we note that we include the so-called 
crystal local-field effects in the normal direction only and neglect them along the surface, i.e., only the 
G = G = element of the noninteracting density response matrix X^tq' (^i ^) has been considered 
in the evaluation of the coefficients Xg=og'=o(^' ^) interacting density response function; 

lateral crystal local-field effects were found to be negligible in the case of the Mg(OOOl) surface ll24l . and 
we expect them to be negligibly small, as well, in the case of other simple metal surfaces like Be(OOOl). 
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Fig. 1 Ab initio band structure of the Be(OOOl) surface. The Fermi level, Ef, is at the zero energy. The grey areas 
show projected bulk Be electronic states. The surface states dispersions are shown by dashed lines. Thin solid line 
shows parabolic fit to the surface state SS located at the F point with effective mass mss = 1.2. The bare surface state 
dispersion with mgg = 1 is shown by thin dotted line. The surface Brillouin zone with two symmetry directions, F M 
and F K, is shown on the right. 

3 Results and discussion 

Figure[T]represents the projected electronic structure of the Be(OOOl) surface, as obtained along the symme- 
try directions F K and T M. The characteristic feature of this surface is the presence of wide energy gaps in 
the projected bulk band structure, where several strongly localized surface states reside ll25l l26l l27l . Here 
we focus our attention on the Shockley surface state that is located around the SBZ center and is marked 
by SS. According to our present first-principles calculations, the binding energy of this surface state at F 
is 2.75 eV in close agreement with the experimental data ll25l l26l l28l |29l [30l and other ab initio calcu- 
lations 1221 EU- As for the energy dispersion of this surface-state band, we find that within the occupied 
part of the band it is almost isotropic and well described by a parabolic-like dispersion with an effective 
mass rrigg = 1.2, again in good agreement with experiment ll26l[32l . However, within the unoccupied part, 
where the surface-state energy approaches the borders of the energy gap, our calculated energy dispersion 
deviates substantially from the quasi-free electron like parabolic behaviour. 

Collective excitations created by an external potential of the form dictated by Eq. (|4]l can be traced to 
the peaks of the absorption probability w(q, uj) of Eq. (|5]l, i.e., to the peaks of the imaginary part of the 
surface response function ^(q, uj) of Eq. (|6]l. We have searched for the maxima of ^(q, oj), as obtained as a 
function of the frequency cu, for various 2D q wave vectors, and we have found the low-energy (acoustic) 
surface-plasmon energy dispersion plotted in Fig.|2]along two distinct symmetry directions: F M and F K. 

As already shown by the ID model potential calculations reported before [11 J, the energy dispersion 
of acoustic surface plasmons follows very closely the upper edge of the continuum (represented in Fig.|2] 
by the grey area) of 2D electron-hole pair excitations occurring within the quasi-2D surface-state band. 
This continuum of electron-hole pair excitations (and, in particular, its upper edge) is, however, dictated 
by the energy dispersion of the surface-state band, which strongly depends on the actual band structure 
and deviates substantially from a quasi-free electron-like parabolic behaviour, as discussed above. This 
results in a first-principles energy dispersion of acoustic surface plasmons that deviates from previous ID 
model potential calculations, shows excellent agreement with the available experimental data along the 
F M direction, and exhibits notable anisotropy. 
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Fig. 2 Calculated acoustic surface plasmon dispersion at the Be(OOOl) surface along two symmetry directions, F M 
and r K, is shown by dots. Grey areas show the "momentum-energy" phase space where intraband electron-hole 
excitations within the SS surface state are available. The experimental data fT2| along the F M direction are shown 
by open diamonds. Thin solid line shows acoustic surface plasmon calculated with the use of the Be(OOOl) model 



potential 1151 and SS effective mass i 



= 1.2. 



We attribute our predicted anisotropy of acoustic surface plasmons to the fact that the upper edge of 
the 2D electron-hole pair excitations occurring within the quasi-2D surface-state band depends strongly 
on the symmetry direction. Our calculations also show that the peaks in the surface loss function are 
stronger and narrower along the F M than in the case of the F K direction. We attribute this fact to the 
complicate Be(OOOl) surface electronic structure, where additionally inter-band transitions involving other 
surface states are also possible (see Fig. [U. It seems that the destructive role of such transitions on the 
properties is more pronounced in the F K direction. Clearly there is no anisotropy in the acoustic surface 
plasmon dispersion obtained with the use of a ID model potential 111], simply due to the isotropic surface 
electronic structure inherent to that approximation. For small momenta, our ab initio acoustic surface 
plasmon dispersion almost coincides with the dispersion obtained with the use of the ID model potential 
and TTigg = 1.2. As the momentum increases, the two calculated curves start to deviate as a consequence of 
the deviation in the corresponding surface state dispersions obtained in the two models. This demonstrates 
the importance of the unoccupied part of the surface state dispersion, as it determines the upper edge for 
the intra-band electron-hole excitation continuum, which represents the crucial factor in the determination 
of the acoustic-surface-plasmon energy dispersion ll33l . 



4 Summary and conclusions 

In summary, we have reported first-principles calculations of acoustic surface plasmons on the (0001) 
surface of Be, as obtained in the random-phase approximation of many-body theory. We have calculated the 
energy dispersion of these collective excitations along the symmetry directions F M and F K, and we have 
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found that the energy dispersion (i) follows closely the upper edge of the continuum of 2D electron-hole 

pair excitations occurring within the quasi-2D Shockley surface state band that is present in this surface, 
(ii) deviates considerably from previous ID model potential calculations, (iii) shows excellent agreement 
with the available experimental data along the F M direction, and (iv) exhibits notable anisotropy. Because 
the energy dispersion of acoustic surface plasmons follows closely the upper edge of the continumm of 
2D electron-hole pair excitations occurring within the quasi-2D Shockley surface state band, we conclude 
that experimental measurements of these plasmons serve to map the anisotropic electron-hole excitation 
spectrum of the surface-state band that is present on the Be(OOOl) surface. The same conclusion should 
be applicable to other surfaces containing a partially Shockley surface state band near the Fermi level, as 
occurs in the case of the (111) surfaces of the noble metals Cu, Ag, and Au. 
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